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We predict the existence of a class of ultracold giant molecules formed from trapped ultracold 
Rydberg atoms and polar molecules. The interaction which leads to the formation of such molecules 
is the anisotropic charge-dipole interaction (a/R 2 ). We show that prominent candidate molecules 
such as KRb and deuterated hydroxyl (OD) should bind to Rydberg rubidium atoms, with energies 
Eh c± 5 — 25 GHz at distances R ~ 0.1 — 1 fim. These molecules form in double wells, mimicking chiral 
molecules, with each well containing a particular dipole orientation. We prepare a set of correlated 
dressed electron-dipole eigenstates which are used in a resonant Raman scheme to coherently control 
the dipole orientation and to create cat-like entangled states of the polar molecule. 

PACS numbers: 33.80.Rv, 31.50.-x , 31.50.Df 



The ability to tailor and control interactions amongst 
species in a quantum gas is a primary motivation for 
topical interface research in atomic, molecular, optical, 
condensed matter and chemical physics. Ultracold polar 
molecules and Rydberg atoms share in this distinction 
for being ideally suited for such studies, as they provide 
a number of convenient control handles [TH5] . Both polar 
molecules and Rydberg atoms can have large permanent 
and induced dipole moments, making them amenable to 
manipulation with external and optical fields via long- 
range van der Waals (R~ 6 ) and anisotropic dipole (R~ 3 ) 
forces. Design of synthetic many-body textures and sim- 
ulations of realistic quantum phases are conditional on 
this level of control [6j [7] . 

It has been proposed that a Rydberg electron inter- 
acting via the zero-range Fermi pseudopotential with a 
nearby perturbing atom could bind into ultralong range 
Rydberg molecular states. Isotropic molecules with MHz 
binding energies, "butterfly" molecules with GHz bind- 
ing energies and large (tens of Debye) permanent dipole 
moments, and "trilobite" molecules with tens of MHz 
binding energies and ultralarge (kDebye) dipole moments 
[H[9] were predicted to exist in ultracold traps. An exper- 
iment with ultracold Rydberg atoms in a dense rubidium 
magneto-optical trap (MOT) recently realized the exis- 
tence of the istotropic utralong range Rydberg molecules 

In this work, we propose the existence of a class of 
ultralong range macroscopic Rydberg molecules which 
can be created in ultracold traps. The molecular bind- 
ing comes courtesy of the interaction of the Rydberg 
electron with the permanent dipole of the trapped po- 
lar molecule. This anisotropic interaction which has the 
form (a — \ )/R 2 is by far, the longest range interaction in 
the physics of cold neutral atoms and molecules and may 
open new vistas into the emergent discipline of long-range 
correlated quantum systems. These giant molecules have 
deep binding energies, ~ 20 GHz, are ultralong range, 
R ~ 2000 ao (see Fig. [2J, and are well protected from 
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FIG. 1: A schematic of the two ultralong range Rydberg 
molecular states which form from the interaction of a Ry- 
dberg atom with a polar molecule. The Rydberg electron 
orbit is exaggerated to highlight the long-range nature of the 
interactions. The direction of the dipole is from the + to — . 



the non-adiabatic avoided crossings, which could lead to 
their premature demise. The double- well structure (left 
(L) and right (R) wells in Fig. in which the orien- 
tation of the polar molecule dipole, is controlled by 
the Rydberg electron, resembles the double-well struc- 
ture which determines the left- and right-handedness in 
a chiral molecule, as shown schematically in Fig. [I] We 
propose a realizable coherent Raman scheme to control 
the orientation of the molecular dipole and prepare cat- 
like superposition of the dipole states. 

The electron-dipole interaction (a — \)/R 2 becomes 
critical at a = 0, resulting from a critical dipole, d cr = 
0.639315 a.u. [11 . For a > 0, the interaction is attrac- 
tive (d > d cr \ and the electron binds to the dipole; for 
a < 0, the dipole is subcritical (d < d cr ) and the electron 
does not bind to the dipole. When d > d cr — 1.63D, 
electron exchange occurs [T2] . 
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For the discussion outlined here, we treat the molecu- 
lar dipole to be subcritical, d < d cr: such that degener- 
ate perturbation theory can be employed to obtain the 
Born-Oppenheimer (BO) potential energy curves. For 
repulsive interactions between the Rydberg electron and 
the molecular dipole (d < d cr ), details of the short range 
physics are unimportant. 

The starting point of our calculation is 




where R n i(r) are the Laguerre polynomials and (....) are 
3 — j symbols, and Z = R Z. 

A molecule becomes "polar", when the two opposite 
parity states of the molecule are mixed to create the per- 
manent dipole moment do. The minimum electric field 
F c required to mix these states with a splitting A is 
F c = Defining the eigenvalues of the potential ma- 
trix as V*(R), the BO potentials for the interaction of a 
Rydberg atom with a polar molecule of dipole do is 



Vx (R) = d \F C - ^ (f* (R) - j^j +F?j (4) 

where F^(R) = V*(R)/d is the electric field at position 
R due to the Rydberg electron and — is the poten- 
tial generated by the ion core. This post-diagonalization 
creates a set of dressed correlated electron-dipole eigen- 
states, with important implications for the control of the 
dipolar superposition. 

The BO potential curves for Rb(n=35)OD( 2 n 3/2 ) 
are shown in Fig. [2] Similar curves exist for 
Rb(n)KRb( 1 H 1 / 2 ) molecule. The modulations are due 
to Rydberg electron oscillations. The interaction of the 
ionic core with the dipole, — , is seen in the potentials 
in which the molecular dipole points toward the core. 
The right panel in Fig. [2] illustrates the transition of the 
molecular dipole orientation through the avoided cross- 



where R is the dipole position with respect to the Ry- 
dberg core and r is the electron position. If we assume 
that the dipole is on the z axis, the matrix elements of 
this interaction in basis \nlm > are 
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ings. The classification of the dipole moments accord- 
ing to their orientation is evident in this figure. The 
dashed curves show the evolution of the dipole which 
points away from the core to that which points toward 
the core (solid curves). The total molecular symmetry 
for the polyatomic molecule is either 2 E or 3 S, depend- 
ing on whether the polar molecule is 1 S, as in KRb, or a 
doublet molecule, as in OD. 

The formation of the giant Rydberg molecular state 
is a competitive process between the electron- and core- 
dipole interactions. Binding occurs when the electronic 
distribution is sufficiently concentrated, allowing the 
electron-dipole interaction to overcome the core contri- 
bution. The concentration of the electron wave function 
is in part controlled by the range of I values available in 
the degenerate Rydberg manifold: I = l m [ n , ...,n — 1. If 
the minimum angular momentum is too large compared 
to n, the electron interaction can not overcome the strong 
1/R 2 core attraction and hence no binding occurs. We 
have found that for l m [ n > 2, n ~ 20 gives an approxi- 
mate cutoff below which no bound states can form. 

The electronic density in the xz-plane at R = 2300 ao is 
shown in Fig. [3] The orientation of the molecular dipole 
is at the top of the figure. The Rydberg oscillations are 
visible with large amplitude at enormous distances from 
the Rydberg core. The presence of a permanent molec- 
ular dipole distorts the spherical symmetry, illustrating 
the anisotropy of the electron-dipole interaction. The 
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With some algebra, the final form for the potential matrix becomes: 
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FIG. 2: The BO potential curves (left panel) for 87 Rb(n=35)- 
OD( 2 II) are shown. The zero of energy is the Rb(n=35) 
threshold. Only those angular momentum states with neg- 
ligible quantum defect (I > 3) in Rb are included. The 
modulations in the potentials are due to the Rydberg oscilla- 
tions and the avoided crossings are due to interactions in the 
degenerate manifold. The two lowest potential wells, near 
R rsj 2000 — 2500 a.u., are weakly coupled with the R (L) well 
supporting vibrational states of the Rb-OD molecule oriented 
away from (toward) the ion core. Positive values of the dipole 
(right panel) refer to the orientation of the molecular dipole 
away from the ion core and the negative values indicate the 
reverse. 

electron wave function seen in Fig. [3] clumps near the 
"positive" side of the dipole with a strong peak in the 
electron probability distribution. 

The size and energy scales of the giant molecules 
generally follow the Rydberg scaling with n. Specifi- 
cally, the outer wells in the lowest two R and L po- 
tentials have minima, respectively, at R e = 1.92 (l)n 2 
and R e = 1.65(2)n 2 with corresponding potential depths 
V R = -0.224(l)<io/n 3 and V L = -0.234(1)<V™ 3 - For 
d < d cr , the potential depth is less than 10% of Rydberg 
energy spacing, allowing us to use degenerate perturba- 
tion theory. At R ei the immense charge separation cre- 
ates a large permanent dipole moment which scales as 
dRm = 1.75n 2 - for n = 35, it translates to dR m = 5.5kD. 

The two outer wells, each corresponding to the molec- 
ular dipole oriented toward or away from the ion core, are 
shown more clearly in Fig. [4j The two wells support a 
large number of vibrational levels within a typical transi- 
tion frequency in the 400 — 500 MHz range. The coupling 
between the L and R curves is on the order ofO(^), 
which for OD is ~ 10 -15 a.u. The diagonal corrections 
to the energies in the n = 35 manifold, are, respectively, 
3 x 10- 9 and 10~ 9 a.u. 

The choice for a realizable polar molecule depends 
on two factors: such molecules should have subcriti- 
cal dipoles, but be of sufficient strength for meaning- 
ful electron-dipole interaction, and have energy split- 
tings which are a fraction of the Rydberg separation. 
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FIG. 3: (Color Online) A density plot of the electronic wave 
function for an n = 35 Rydberg molecule at R = 2300ao is 
shown in cylindrical coordinates. Darker areas show larger 
amplitudes of the electronic wave function. The position 
of the Rb + (red circle) and the orientation of the molecular 
dipole (green-blue dots) are indicated. The size of the diple is 
greatly exaggerated for illustration. The spherical symmetry 
is broken by the anisotropic electron-dipole interaction. This 
density plot is a realization of the right configuration in Fig. 

m 

The hydroxyl radical (OH), a A-doublet molecule ( 2 n 3 / 2 , 
Aa = 1.67 GHz), was recently magnetically trapped [T3] 
and its cold collision with Xe and He beams were studies 
[14]. Deuteration of hydrogen radicals reduces the dou- 
blet splitting by several factors [16], making them more 
readily polarizable. 

In Table [I] we give a list of favorable polar molecules. 
The ground state of 40 K 87 Rb has been produced with 
large phase space density and its dipole moment and ro- 
tational splitting {B e ) have been established [2]. While 
the dipole moment for the KRb molecule is small, its 
small rotational splitting {2B e ) : and its ubiquity make it 
a good candidate. The deuterated radical, CD( 2 n 3 / 2 ), 
can be polarized with field strength of F c ~ 2 x 10 -7 
a.u. and has a large but subcritical dipole moment, 
1.46 D [15]. Another potentially realizable candidate is 
metastable CO(a 3 ni 52 (^ = 0)) which lives for 2.63 ms 
and its A± transition (O = 1, J = 1) is 394.1 GHz and 
has a dipole moment of 1.38 D [T7] . 

The lowest vibrational levels in L (0l) and R {Or) 
wells in Fig. [4] can be coherently coupled in a microwave 
Raman process. This fully on-resonance, but weak in- 
tensity two-photon transition, via an intermediate vi- 
brational state {vr) in the extended right well (dashed 
curve) allows for efficient transfer of population. The 
configurations in the L and R wells, respectively, are with 
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TABLE I: Choice of polar molecules and their properties. 



1 cm 



: 30 GHz. 



molecule 


d (D) A(MHz) 


F c (a.u.) 


B e (cm -1 ) 


Ref. 


od ( 2 n 3/2 ) 


1.66 494 


6 x 10 _B 


-40 


E2 


KRb (%) 

cd ( 2 n 1/2 ) 
CO (a 3 n 1)2 ) 


0.566 


6 x 10~ 7 


0.0371 


m 


1.46 897 


1.6 x 10" 7 


- 60 


15 


1.38 394 


7 x 10~ 8 


40 


EH 




ization fields, on-resonance microwave field strengths of 
Ffj, < F c are desired. With A± = 20, where Q = d tr • 
is the Rabi frequency, microwave fields of ~ 10 -7 
a.u. (500 V/cm) can simultaneously maintain coherence 
in the cat states and drive the Raman transitions. 

In summary, we explore the possibility that giant poly- 
atomic Rydberg molecules could readily form in ultra- 
cold trapped mixture of atoms and molecules, from in- 
teraction of Rydberg atoms and polar molecules. These 
molecules which form in double wells contain dipole con- 
figurations which can be coherently controlled by a Ra- 
man process. The dipolar interaction of two polarized 
giant Rydberg molecule could be coupled to the internal 
cat states of the polar molecules. Ionization of the Ryd- 
berg atom from each L and R well can deterministically 
measure the orientational state of the polar molecule. 

We are grateful to M. Cavagnero, C. Ticknor and J. 
Feist for discussions. This work was funded by a grant 
from the NSF to ITAMP at the Harvard- Smithsonian 
Center for Astrophysics and Harvard Physics Depart- 
ment. 



FIG. 4: The two (L and R) wells of the Rb-OD Rydberg 
molecule are shown. Also shown, are a select number of vi- 
brational energy levels and associated wave functions, and 
our proposed Raman scheme for coherent population trans- 
fer between the two wells and preparation of cat-like states. 
The quantum numbers refer to the lowest vibrational states 
in the left (Oz,), right (Or), and an excited level in the right 
well (vr), and Ql, and Qr are the Rabi frequencies for the 
Raman transitions. 

the molecular dipole pointing toward or away from the 
ionic core. The coherent Raman scheme ensures that the 
molecular dipole orientations become entangled. Since 
the distance between the dipoles scales roughly as 0.27n 2 , 
which for n = 35 translates to 16.5 nm, cat-like superpo- 
sition states are possible. 

We compute the dipole transition strengths for the Ra- 
man process. Due to the existence of an inner wall in the 
potential curve which contains the R well, large Franck- 
Condon (FC) factors between the L and R states with 
the upper vr Raman states can be found. FC factors for 
0l —> vr and Or — » vr transitions are typically about 
0.02 a.u. The electronic transition dipole moments, with 
photons polarized along the internuclear axis, scale as 
0.20n 2 ; for n=35 Rydberg state, the transition dipole 
has a value of 265.6 a.u., giving a total transition dipole, 
d tr — 5 a.u. To preserve the coherence of the cat states, 
\c >= -jz[\R > -\-e l ^\L >], in the presence of the polar- 
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